Entropy scaling laws for diffusion 
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Recently, Samanta et al. Q set out to understand the low density failure of the entropy scaling law for the self- 
diffusion coefficient D conjectured by Dzugutov and to provide a simple alternative. After an interesting derivation, 
, that however contains a number of uncontrolled approximations, they arrive at Eq. 7 of 0, which reduces for a hard 
sphere fluid to: 



" ' (1) 



De 1-Se/fcs 

O with De = Db/x ^^^d Db = 3{kBT/Trm)i /8pa'^ the Enskog and Boltzmann diffusion coefficients "s], x the contact 
^ value of the pair correlation function, Sg the excess entropy per particle, and A = 2.5. Using the excess entropy 
^ ■ Se/kB — — (4f7 — 3?7^)/(l — 77)^ and contact value x = (2 — r])/2{l — rj)^ given by the Carnahan-Starling equation of 
^l* I state (t? = Trpa^/6 - packing fraction), we test this relation against the molecular dynamics simulation results of 
^ • Erpenbeck and Wood (EW) 5] - Fig. 1. The disagreement is quite significant and more important the behavior of the 
_ two curves is very different, as {D/De)ew ^^"-"^ ^ monotonically decreasing function of {^s^/kB). The discrepancy 
cannot be attributed to the authors use of an approximation for the excess entropy. 

The idea of entropy scaling for transport coefficients has a fairly long history 0. Arguing on the basis of the 
molecular "caging" effect Dzugutov (2J proposed the scaling law: 



g: = Bexp(se//cB) (2) 



which assumes that the natural length and time scales for diffusion are given by a suitably defined hard sphere diameter 
cr and the Enskog collision frequency Te — 4:a^xP\/^kBT/m, with B an universal constant. Unfortunately, the above 
relation appears to work only in a limited, high density domain for both hard spheres j?! and realistically modeled 
fluids 8J. The problem that arises at low and moderate densities (see Fig. 1) with the scaling introduced in Eq. 2 can 
10 , be understood if we observe that the left-hand-side of that equation can be written up to a multiplicative constant 
CO I as D / DbXtJ^ ■ Therefore, in the limit of a dilute system, 77 ^ 0, this term diverges as l/rj^, while the right-hand-side 
of Eq. 2 remains finite. This behavior should be expected for any valid definition of a and x Se approximation. 
The noted pathology of Eq. 2 can be avoided by replacing a as the relevant length scale with l/pcr^, the Boltzmann 
' mean-free path, which should be a reasonable measure of the degree of molecular confinement, oc a/rj. The 

new relationship is: 



^ exp (jSe/kB) (3) 



DbX 

^ . . 

(~| where we introduced a different constant 7. The test of this dependence is shown in Fig. 1 for hard spheres, with 

O ■ 7 = 0.8. Furthermore, Eq. 3 holds for Van der Waals fluids as well f^. 

^ ' Samanta et al. also propose a generalized Stokes-Einstein relation. However, such a relation is hardly necessar 
^ I given that the usual Stokes-Einstein formula with the 'slip' boundary condition holds well for both hard spheres 
and Van der Waals fluids . 

This work was performed under the auspices of the U. S. Department of Energy by University of California Lawrence 
' Livermore National Laboratory under Contract No. W-7405-Eng-48. 
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FIG. 1: Comparison of the hard spheres diffusion coefficient (Ref. Q) - circles, with scaling relation of Ref. Ifj (Eq. f) - dashed 
line, Dzugutov scaling law |2| (Eq. 2) - dot-dashed line, and new entropy scaling (Eq. 3) - solid line. 
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